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Abstract 

A system with anholonomic constraints where the trajectories of physical degrees 
of freedom are autoparallels on a manifold equipped with a general Cartan connec- 
tion is discussed. A variational principle for the autoparallel trajectories is derived 
from the d Alambert-Lagrange principle for anholonomic constrained systems. A ge- 
ometrical (coordinate-independent) formulation of the variational principle is given. 
Its relation to Sedov's anholonomic variational principle for dissipative systems and 
to Poincare's variational principle in anholonomic reference frames is established. A 
modification of Noether's theorem due to the torsion force is studied. A non-local 
action whose extrema contain the autoparallels is proposed. The action can be 
made local by adding auxiliary degrees of freedom coupled to the original variables 
in a special way. 

1 Anholonomic constrained systems 

There is no need to explain how important are constrained systems in modern physics 
(e.g., electrodynamics, Yang-Mills theory, general relativity, etc). Constraints in dynami- 
cal systems are usually regarded as a part of the Euler-Lagrange equations of motion which 
do not involve time derivatives of order higher than one. In other words, both constraints 
and equations of motion result from the least action principle applied to some Lagrangian. 
The existence of the Lagrangian formalism is of great importance in constrained systems 
because it allows one to develop the corresponding Hamiltonian formalism 0] and canon- 
ically quantize the system [£]]. Yet, the variational principle is a powerful technical tool 
to find integrals of motion of dynamical systems via symmetries of the Lagrangian. 

The Hamiltonian (or Lagrangian) constrained systems form a relatively small class 
of constrained dynamical systems. Given an "unconstrained" system whose dynamics is 
governed by a Lagrangian L = L(v,x), v l and x % being generalized velocities and coordi- 
nates, respectively, one can turn it into a constrained system by imposing supplementary 
conditions F a (v,x) = (constraints) which has to be fulfilled by the actual motion of 
the system. There two ways to incorporate the constraints into a dynamical description. 
First, one can simply modify the Lagrangian L —> L + X a F a with X a being the Lagrange 
multipliers and treat the latter as independent dynamical variables. In doing so, we are 
led to the Lagrangian constrained dynamics. The other way is to supplement the uncon- 
strained Euler-Lagrange equations d/dt(d v L) — d x L = by the constraints F a = 0. It 
is well known that if the constraints are not integrable, the two dynamical descriptions 
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are not equivalent [|2], H. The non-integrable constraints are called anholonomic, and the 
dynamical systems described in the latter of the above two ways are known as anholo- 
nomic systems. There exists no Lagrangian or Hamiltonian formalism for anholonomic 
systems, i.e. they are non- Lagrangian dynamical systems ||. The existence of constraints 
implies that the dynamical system has non-physical degrees of freedom, meaning that the 
actual motion of the system is determined by a less number of independent parameters 
than the number of initial values of generalized coordinates and velocities. For holonomic 
systems, the motion of the physical degrees of freedom can be obtained by applying the 
conventional variational principle to the Lagrangian reduced on the constraint surface 
L — > L\f = q. When applied to anholonomic constraints, this procedure leads to wrong 
equations of motion 0, However dynamics of physical degrees of freedom in anholo- 
nomic constrained systems may also possess "good" properties, e.g., a covariance under 
some group transformations, existence of integrals of motion with a clear physical inter- 
pretation, etc, that is, the properties that one always wants to see in physical systems. 

An example of this kind is provided by the autoparallel and and geodesic motions on a 
manifold equipped with a general connection T^ va compatible with metric g^ u , D^g ua = 0, 
where is the covariant derivative. In it has been shown that the autoparallels can 
be realized as the trajectories of the physical degrees of freedom in a special anholonomic 
constrained system, while the geodesies can always be regarded as the trajectories in a 
holonomic constrained system. On a manifold M the autoparallels and geodesies are 
determined respectively by the following equations covariant under general coordinate 
transformations 

D v v^ = ti^ + rV/^O; (1.1) 
D v v» = v» + r» ua v v v a = . (1.2) 

Here D v is the covariant derivative along the velocity vector v M and T fJi ucr are the Christoffel 
symbols. Both Eqs. (|1 . 1| ) and (|1.2|) determine a curve that parallel-transports its tangent 
vector along itself. The curve with such a property is the autoparallel when the most 
general connection compatible with the metric is used to specify the parallel transport. 
The geodesic motion occurs if the natural Riemannian connection (induced by the metric) 
is chosen to define the parallel transport. The difference of ( |L1| ) from ( |1.2| ) resides in the 
torsion force term. Any connection compatible with the metric can always be represented 
in the form [|5J gaX^^^u = I" apiv 

+ K a ^ u , where K aiiV = S afiu - S^ a + S vrTtx is called the 
contorsion tensor, and g aX S\ fJjU = S a ^ u = |(r CT M „ — T a Ufl ) is the torsion tensor. The 
deviation of the autoparallels from the geodesies is caused by the torsion force K^ va v v v u . 

The geodesic equation ( |1.2| ) follows from the Hamilton variational principle applied to 
the action 

S=~Jdtg^v 1/ , (1.3) 

or to its relativistic analog 

S= - I ds = - f dty/g^WV , (1.4) 
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with s being the proper time (or length) of the trajectory (in this case, t> M = dq^/ds and 
■u M = dv^/ds in ( |1.2|) ). A particle moving along the autoparallel trajectory is an example 
of a non-Lagrangian system. A system of differential equations of second order is called 
Lagrangian if there exists a local Lagrangian such that the original system is equivalent 
to the Euler-Lagrange equations. One can formulate the sufficient conditions for a given 
dynamical system to be Lagrangian || . Even if these conditions are not fulfilled, one can 
still try to find a non-local action for a given dynamical system. In Section 6 we construct 
an explicit example of a non-local action whose extrema contain the autoparallels. From 
the geometrical point of view Eq. Q) is just as good as Eq. (|1.2|) and may be regarded 
or postulated as an equation of motion of a spinless particle on a manifold Q. 

Here we generalize the embedding procedure of [|J to arbitrary spaces with curvature 
and torsion. Next we make use of the representation of the autoparallel motion as a 
motion of an anholonomic system to establish a variational principle for the Lagrangians 
( |1.3j ) and ( YA ) which leads to Eq. fll.lp . The variational principle is derived from the well- 
known variational principles for anholonomic systems, such as Gauss' principle of least 
constraint, Holder and d'Alambert-Lagrange principles [FJ, ^|, |TU|. But in contrast to them, 
it has an advantage that it applies to Lagrangians reduced on the surface of constraints. 
We shall also show that the new variational principle can be given a completely covariant 
(coordinate-independent) formulation on a manifold with a general Cartan connection. 
For this reason we shall refer to it as a covariant variational principle. Its relation to 
the variational principle of Poincare [0] and to that proposed by Sedov for dissipative 



systems |12j is explained. Finally, we propose a modification of the actions (|1.3| ) and ( |L4|) 
by adding new auxiliary degrees of freedom so that the modified actions have extrema 
being the autoparallels and admit the conventional Hamiltonian formalism. 



2 Autoparallels from constrained motion 

Consider a metric manifold M and local coordinates g M on it. Let T^ ua be components of 
a connection on M in the coordinate basis. We denote V(M) the space of all paths in M, 
T q M the tangent space at a point and TM the tangent bundle. Consider an auxiliary 
Euclidean space H n of the dimension greater than that of M, n > dimM. Cartesian 
coordinates in IR n are denoted by x % . In the space "P(IR n ) of all paths in H n we define a 
subspace of conceivable paths (i.e., of those allowed by constraints) as an image of V{M) 
in the embedding V(M) -> V(JR n ): 

**(«) = fdq^iq) , (2.1) 

for any path g M (s) in M. The embedding functions e l (q) are smooth on M. From the 
path embedding ( |2.1| ) follows the embedding of the tangent space T q M into IR™ 

v l = 4(gK , v»eT q M, (2.2) 

2 Here we do not discuss physical arguments which of those two trajectories should be identified with 
a physical trajectory of a spinless point-like particle moving in a space with torsion. An interested reader 
may find physical arguments supporting both of them, geodesies Q and autoparallels J|, [J . 
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because one can always find a curve passing through a point such that t> M = 

dq»(s)/ds. 

The space M can not be embedded into IR n pointwise if the constraints ( |2.2|) on 
the tangent space (or on the velocities of the conceivable motion) are not integrable 
(anholonomic constraints): 

<fdq%(q)^Q, (2.3) 

for any closed path in M. That is, the one-forms e^dq^ are not closed, and there exist no 
mapping M — > IR™, x % = x l {q), induced by ( |2.2|) . 

The scalar product in IR™, (w,v) = 5ijW % v^ induces the metric on M. For any two 
vectors being the images of two elements of T q M we have 

(w, v) = g^wW , = (e„, e v ) , v", w" G T q M . (2.4) 

Next, we determine a connection induced by the path embedding. To do so, we require 
that the parallel transport in M must be compatible with the embedding ( |2.1|) . This 
means the following. Consider a curve q fl {s) passing through a point g" G M and its 
image x l (s). Take a vector field G T q M and parallel-transport it along the curve 
g M (s). The resulting vector is then embedded into IR n by means of (|2.2|) . Let us repeat 
the procedure in the opposite way. We first embed the vector u> M and then parallel- 
transport its image along x l (s), the image of g M (s). The compatibility condition implies 
that the vectors obtained by these procedures may differ by a vector orthogonal to the 
hyperplane being the image of T q M, and this should hold for any curve and any w^. We 
remark that this condition is weaker than a similar condition of |3| and it is sufficient to 
obtain the most general spaces with curvature and torsion. 

An infinitesimal variation of under the parallel transport is proportional to D v w^, 
where t> M = g M , while the corresponding variation of its image, w l = e^w 11 , is D v w l = 
d v w % = dw l /ds because the connection in IR™ is trivial. The compatibility condition 
means that 

dm 1 

— - e\D v w^ = f; u vW , (e a , f, v ) = . (2.5) 

Computing the derivative in the left-hand side of (|2.5|) and assuming v M and w^ 1 to be 
arbitrary we find 

D,ei = . (2.6) 
Multiplying this equation by e % a we obtain the connection coefficients 

r^ = ^ A (£ A) cU CT ) , (2.7) 

where g v ^g IMJ = 5%. It is easy to verify that the compatibility condition holds for tensors 
of higher rank, too, provided the induced metric is used to lower and rise tensor indices. 
Thanks to the compatibility condition, the induced connection turns out to be compatible 
with the induced metric D^g uu = (D^e u ,e a ) + (e„, D^e a ) = 0. In the non-integrable case 
( |2.3| ), the induced connection may have a torsion 

S\a = \f X [(e x , d v e a ) - {ex, d a e v )\ . (2.8) 
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The tensor fields /* must also satisfy the integrability condition for Eq. ( |2.6| ) resulting 
from the commutation relation 

[Dfi, D v ] = R^ aX L aX - 2S% lu D a , (2.9) 

where = —L Vil are generators of a (pseudo)orthogonal group acting in the tangent 
space of M, and R flu crX is the Riemann-Cartan curvature tensor. Making use of the 
integrability conditions we derive 

{ffiX, fucr) — (fuX, f[ia) = ^R^uXa ■ (2-10) 

That is, the tensor fields j % determine the Riemann-Cartan curvature tensor on M. 
The case f 1 = corresponds to the so called teleparallel spaces with zero curvature and 
nonzero torsion. With sufficiently large n (at least one should have In > l + (dimM) 2 , the 
number of independent components of the metric and torsion tensors), one can construct 
the most general connection on M. 

Let us analyze dynamics induced by the anholonomic constraints ( |2.2|) , assuming the 
unconstrained motion to be a free motion in the Euclidean space H n . There are several 



variational principles for anholonomic systems [1C] that can be applied to our system to 
derive the equations of motion. If the constraints are integrable, they all are equivalent 
to the conventional Hamilton variational principle. 

In the auxiliary Euclidean space the states of the system is labeled by pairs = (v \ x l ). 
For any two trajectories x\ 2 (s) passing through the state ip, Gauss' deviation function is 
defined as || | 

G^ = \(v\-vi)H t3 (v{-vi) , (2.11) 

where Hij(ip) = d 2 L/dv l dv^ is the Hessian for a Lagrangian L in the state ip. Gauss' 
principle of least constraint says that the deviation of conceivable motions (allowed by 
the constraints) from the released (unconstrained) motion takes a stationary value on the 
actual motion. A physical meaning of Gauss' principle is transparent: the acceleration (or 
the force) caused by the constraints must have a minimal deviation from the acceleration 
of the unconstrained motion. In our system the released motion is the free motion v l = 
0,Hij = the accelerations of the conceivable motion are obtained by taking a time 
derivative of ( |2.2j ). The Gauss deviation function assumes the form 

GV = \g^D v v v D v v^ + ^(f v , f v ) , /* = f^v" . (2.12) 

The second term in the deviation function does not depend on the acceleration at the 
physical state xjj = (g^t^), while the first term is non-negative and attains its absolute 
minimum if D v v^ = 0, i.e., for the autoparallel trajectories. 

Let us denote the Lagrange derivative as [L]i = —d/ds(d v tL) + d x iL. We remind that 
elements of the tangent space are also called virtual velocities in analytical mechanics. 
The d'Alambert-Lagrange principle asserts that the conceivable motion of a system with 
the Lagrangian L is an actual motion if for every moment of time ||. 

(w,[L]) = (2.13) 
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for all virtual velocities of the constrained motion. We take L = (v,v)/2 and calculate 
[L]i for the conceivable motion [L]j = — v % = —d/ds{e l ^v^). For virtual velocities of the 
constrained motion w l = s^w^, Eq. ( |2.13| ) assumes the form (w, [L]) = —g liu w u D v v^ = 0. 
For an arbitrary w 11 , it leads to the autoparallel equation Ql.lj) . 

There is an equivalent formulation of the d'Alambert-Lagrange principle known as 
Holder's variational principle [§, 0]- A conceivable path is called a critical point of the 
action functional if its variation vanishes when restricted on the subspace of virtual ve- 
locities of the constrained motion. Holder's variational principle assumes that the actual 
constrained motion is the critical point of the action. For an infinitesimal variation of the 
trajectory we have 5x z (s) = uw l (s),u — > 0. Then SS[x]/Sx % = d u S[x + 5x]\ u =o = 0. This 
is equivalent to (|2.13|) . Restricting the virtual velocity w l to the subset specified by the 
constraint ( |2.2| ) we arrive at the autoparallel equation ( |1 . 1| ) . 

All the conventional variational principles for anholonomic systems are applied to non- 
constrained Lagrangians, while the constraints are implemented through the restriction of 
path variations to a specific class determined by the constraints. We shall now develop this 
idea further and find an equivalent variational principle that applies to the Lagrangian 
restricted on the surface of constraints. We recall that for anholonomic constraints the 
Hamilton variational principle leads to wrong equations of motion, when applied to the 
Lagrangian restricted on the constraint surface. In particular, we shall find a variational 
principle for the Lagrangian (|1.3|) or (|1.4|) that leads to the autoparallel equation (|1.1| ). 

Before we proceed let us make a remark. Our treatment of the autoparallel motion as 
an anholonomically constrained motion is somewhat unconventional. Typically, anholo- 
nomic constraints are imposed directly on positions and velocities of the unconstrained 
motion, i.e. on x % and v l (e.g., a skater on an inclined plane, a rolling ball on a rough 
surface, etc, [|3], |T3]]). They can be regarded as restrictions on the initial values of veloci- 
ties and positions so that solutions of equations of motion depend on the less number of 
parameters. A historical reason for such a treatment is that the positions x l can be mea- 
sured and have a natural physical meaning in many concrete anholonomic systems [13| . It 



is important to observe that for anholonomic constraints there exists no subspace in the 
original configuration space that could be identified with the physical configuration space 
formed by initial values of positions allowed by the constraints. In our anholonomic sys- 
tem, we have explicitly introduced a physical configuration space M, while the auxiliary 
Euclidean space is used only to formulate a dynamical principle that specifies the motion 
in M. Since anholonomic constraints do not allow us to regard M as a submanifold in 
the configuration space IR n of the unconstrained system, we have imposed constraints 
by restricting paths in R" to a subclass permitted by the constraints, i.e., through the 
embedding of all trajectories of physical degrees of freedom V(M) into V(JR n ). This is 
always possible for any type of constraints. 

3 Variational principle on manifolds 

Here we shall give a formulation of the variational principle on manifolds that is convenient 
for the subsequent generalization to anholonomic systems. 
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Consider a vector field w^{q) on M satisfying the boundary condition 

w »( qi ) = w ^q 2 ) = , (3.1) 

and an action functional S = J dsL(v, q) for any trajectory connecting the points q± 2 . We 
define a variation of the action relative the vector field by 

5 W S = Jds i^—d^ + — w " j . (3.2) 

The derivative d w v M of the velocity along the variation vector field specifies the variation 
of iA Given a trajectory q^ = g M (s), the velocity vector field v M is known only along the 
trajectory, while we need to know the behavior of v M in the vicinity of the trajectory in 
order to calculate d w v M for a generic We require 

C w v» = d w v^ - d v w^ = , (3.3) 

where C w v^ = —C v w^ denotes the Lie derivative on M. Thus, for any Lagrangian, being 
a function on the tangent bundle TM, the smooth vector field w^ 1 on M determines a 
variation of the position, and Eq. ( |3.3|) specifies the variation of the second independent 
coordinate on TM, the velocity v fl . Eq. ( p.3| ) can be given a tensor form symmetrical 
relative the velocity and variation vector fields 

CvV* = C v w» . (3.4) 

The Euler-Lagrange equation follows from d v w^ = dw^/ds, the boundary condition ( |3 . 1|) 
and S W S = that should hold for any w^. Since the Lie derivative of a tensor is a tensor 
0, Eqs. ( |3.3| ) and ( p.4[ ) are covariant under general coordinate transformations, and so 
are the corresponding equations of motion if the Lagrangian is a scalar. 

The above geometrical formulation of the variational principle is equivalent to the 
conventional one where the variation of the action is defined via smooth path variations. 
Indeed, equation ( |3.3| ) implies that for any trajectory g M (s) with fixed endpoints there 
exists a one-parameter family of trajectories q^ l {s 1 u) = q£(s) with the same endpoints 
such that g^(s,M = 0) = q^(s) and = d u q^(s,u) in the vicinity of the trajectory g M (s). 
In other words, there exists a local coordinate net g M (s, u) such that v M and w 11 are tangent 
vectors for the coordinate lines u = const and s = const, respectively. The variable u 
plays the role of the variation parameter so that 5 w S[q) = d/duS[q u ]\ u =o. 

If the velocity variation is specified by means of a new principle, other than a smooth 
deformation of the path, the condition ( |3~3| ) would get modified and, therefore, the vari- 
ational principle would yield new equations of motion. To obtain such a new principle 
for the autoparallel motion, we make use of the variational principles for anholonomic 
systems discussed in Section 2. If the condition (|3.3| ) is dropped, then there exist no 
coordinate net q^(s,u) such that u> M = duq^i^s^u) if = dsq^^s, u). For this reason the 
variation ( |3.2|) will be called non-coordinate or anholonomic. Note that remains a 
smooth vector field on M. 
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4 Covariant variational principle 

On a path with fixed endpoints, consider a vector field u> M that satisfies the boundary 
condition (|3.1| ). Its image w l (s) = determines a variation vector field on 



the image trajectory x l (s). By construction the variation vector field w l (s) belongs to the 
class of virtual velocities allowed by the constraints. The variation of the velocity vector 
d w v^ is to be determined by Holder's variational principle ( |2.13| ). We assume that the 
variation (|3.2j ) of the unconstrained action (fi is to be replaced by the Cartesian index i) 



is equal to the Holder variation ( |2.13| ) when the former is restricted on the constrained 
surface 

5 W S= / ds(w,[L]) . (4.1) 



Integrating by parts in the right-hand side of ( |4.1| ) and making use of the boundary 
condition ( |3~TD we obtain 

The integration has been omitted because (|4.1| ) holds for an arbitrary variation vector 
field allowed by the constraints. On the constraint surface, we have d v iL = efd v ^L, = 
Sijg^el and w l = e^w^. Relation (fi~2"l) leads to 

(d w v - d v w, e") = , (4.3) 

since d v u,L is also arbitrary. From the compatibility condition (|2.5| ) follows that 

(d w v,e> 1 )=D w v> 1 , (4.4) 

and, similarly, 

(d v w, e») = D v w>* . (4.5) 
Eq. ( |4.3|) leads to the sought condition that specifies the variation of the velocity v M 

D w v" = D v w" . (4.6) 



It can be written as 



£ wV » = d w v» - d v w^ = 2S^ (Jv v v w' J . (4.7) 



The derivative dwV^ is proportional to the difference of the vector field v 1 * at two neighbor- 
ing points g M and q' 1 +uw fl , u — > 0, i.e., ud w v^ = (v fJ- +ud w v t1 )— v^. A similar interpretation 
holds for d v w^. The left-hand side of Eq. (|4.7| ) contains four vectors that can be combined 
to form a parallelogram [w^ + (v M + ud w v^)\ — [v 1 * + (w^ + ud v w^)] which is not closed 
as follows from (|4.7Q . Thus, Holder's variational principle has led us to the conclusion 
that the velocity variation must be chosen so that the closure failure of the parallelogram 
formed by the velocity and variation vectors would be proportional to the torsion. Note 
also that the closure failure of the parallelogram induced by the parallel transport of any 
two vector fields along one another (|4.6j ) is also used to define the torsion on a manifold 

0. 
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Let us take a Lagrangian L = L(v, q) on TM and find the equation of motion resulting 
from the new variational principle. Substituting d w v M into ( |3.2|) and making use of the 



boundary condition (O) to integrate by parts, we get 

5 W S = J dsw^ ([LI + 2SV|^>) • (4i 



The modified Euler-Lagrange equations are 

Eq. ( |4.9|) has also been considered in [14], [15| to describe the motion in spaces with torsion. 



Taking the Lagrangian L = (v, v)/2 or L = — y (v, v) in the auxiliary Euclidean space 
and restricting it on the constraint surface v l = s^v^, we obtain the action ( |1.3|) or 
( |1.4| ). By construction the variational principle 5 W S = 0, where the variation on TM is 
determined by ( |4.6|) , should yield the autoparallel equation ( |1.1|) . It is not hard to be 
convinced that the modified Euler-Lagrange equation ( (4.9| ) indeed leads to (|1.1|). 

A few remarks are in order. Our derivation of the condition ( f4.6|) does not rely on 
that whether the unconstrained Lagrangian explicitly depends on the auxiliary Cartesian 
variables x l or not. The terms involving d x tL in Eq. ( |4.1| ) are cancelled. For this reason 
the final condition (}4.6|) does not depend on the form of the auxiliary unconstrained 
Lagrangian and may be applied to any Lagrangian on the physical configuration space. 

The condition (|4.6|) is covariant under general coordinate transformations on M and 
has a transparent geometrical meaning. A covariant variation of along is equal to a 
covariant variation of w M along v^. All the modification of (|3.4| ) we have made is that the 
Lie variation has been replaced by the covariant variation. It is quite remarkable that for 
tensors on a manifold there exists only two independent variations that produce tensors out 
of tensors and involve a displacement ||, p. 336: The Lie variation £ w and the covariant 
variation D w . Thus, two geometrically distinguished curves on a manifold, geodesies and 
autoparallels, can be associated with the two independent variations available on the 
manifold equipped with a connection compatible with the metric. 

If the torsion is not zero, the variations induced by operators D w and D v are non- 
coordinate (or anholonomic) because the basis in T q M is a non-coordinate basis fT6]| . 
Indeed, assume that there exists a coordinate net g M = q^(s,u) such that the relation 
d S)U = D v<w holds. Taking F to be a scalar, from (|2.9|) we derive 



[D fl ,D u ]F=-2S% v D a F . (4.10) 



The curvature term does not contribute to the commutator ( 4.10| ) because F is a scalar 



Thus, [d S} d u ] 0, and there is no coordinate system for which the covariant derivatives 
D V)W play the role of the translation operators along the coordinate lines. 

The use of a non-coordinate basis in the tangent space to determine a variation of 
the action is not something unusual in mechanics. Some analogy can be made with 
Poincare's variational principle in non-inertial reference frames (e.g. a rigid body in 
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the body-fixed reference frame) [TT], |3|]. One looks for a re-formulation of the Hamilton 
variational principle in a non- coordinate basis in the tangent space of M. The coordinate 
basis in T q M can always be chosen as <9 M so that a variation of F(q) is ~ d^F. We can 
also assume another basis e M / = e t t,{q)d tl in the tangent space to determine a variation of 
any quantity on M. In general, this basis is non- commutative and, hence, non-coordinate 

[e^,e vl ] = 2C r7 '^ u ,e a , , (4.11) 



where the structure functions of the Lie algebra (|4.11|) are coefficients of the object of 



anholonomity |J. In the new basis we have t> M = e^,f M , e^,e% = 5%. The components 
v M ' of the velocity vector field in a non-coordinate basis are called quasivelocities because 
there is no g M (q) such that = cf . The problem is to find a variational principle 
for the Lagrangian where the velocity components are taken in the non-coordinate basis, 
L{v^ , q^) = Liv 1 *, q^). It is solved by rewriting Eqs. ( |3.2|) and (|3.3| ) in the non-coordinate 
basis. Eq. ( |3.3| ) assumes the form ( |4.6| ) where /i is replaced by // and the covariant 
derivative is taken relative to the connection, r°" M v = e a a e^d^e",, induced by going 
over to the new basis. Such a connection has an antisymmetric part equal to the object 
of anholonomity. The variational principle SS/Sw^ = yields celebrated Poincare's 
equations. They have the form ( }4.9| ) where L — > L, d v ^ — > d v ^ , d ql + — > e^,d q v and S a ^ iu — > 



C u fj,/ u '. 

As we see Poincare's variational principle is based on a coordinate variation (a smooth 
deformation of the path), but the variation of the velocity components is taken in a non- 
coordinate basis in the tangent space. Thus, no torsion force can be gained by considering 
a variational principle in an anholonomically transformed basis. By definition ||, the 
torsion transforms as a tensor S a = e° e^,e^ l S a ^ v = (T a — T a u /^/)/2 — C a 
If the torsion is zero in one basis it is zero in any other. The covariant variational 
principle always induces the torsion force because the condition ( |4.6| ) is covariant under 



all (coordinate and non-coordinate) transformations of the basis in the tangent space. 
But it should be kept in mind that the variation specified by the condition ( fO| ) is no 
smooth deformation of the path with fixed ends. 

As a final remark, we shall point out that anholonomic variations in analytical mechan- 
ics have also been introduced by Sedov |L2| to study dynamics of dissipative systems (they 
are examples of non-Lagrangian systems). He also proposed an anholonomic variational 
principle for such mechanical systems and considered its applications. 



5 Noether's theorem 

In Lagrangian mechanics first integrals of motion can be obtained from Noether's theorem. 
The covariant variational principle has led us to the new equations of motion ( [4 . 9| ) . The 
presence of the torsion force should affect Noether's integrals of motion. Therefore it is 
natural to expect Noether's theorem to be modified. 

Consider a one-parameter group of diffeomorphism on a manifold. Given a trajectory 
g M (s), a vector field determines a smooth deformation of the trajectory under the 
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one-parameter group of diffeomorphism on the manifold 



d UJ q» = u 1 *, d u) v^ = u' t . (5.1) 

Let the Lagrangian be invariant under the transformations ( |5.1| ) up to a total time deriva- 
tive 

d u L = —, $ = $(v,q,s) . (5.2) 
as 

If the motion is determined by the Euler-Lagrange equation, then in accordance with 
Noether's theorem ||, the system possesses the integral of motion 

§ = »■ ' = ^ 

The proof follows from relation (|5.2|) that should be written in the form 

f g + [L],^ = . (5.4) 

For the actual motion [L] M = 0, thus leading to the conservation law (|5.3|) . Similarly, 
expressing [L]^ from ( f4.9| ) and substituting it into ( |5.4j ), we derive modified Noether's 
theorem 

fjUMvH^. (5.5) 
In particular, for the time translation symmetry we have lo^ = and $ = L. The right- 



hand side of (5J3) vanishes since the torsion tensor is antisymmetric. The corresponding 
integral of motion is the system energy. On the other hand, assuming the Lagrangian to 
be invariant under spatial translations and rotations (e.g., L = |t> 2 , g^ v = 5^ v ), we observe 
that the momentum and the angular momentum are no longer integrals of motion for a 
generic torsion tensor. 

One can give two additional equivalent formulations of Noether's theorem. Assume 
that under the transformations ( |5.1| ) the following relation holds 



Then ( |5.3| ) is an integral of motion. Clearly, to achieve ( |5.6|) , the vector field should, in 
general, depend on the torsion. Although the torsion force violates the Noether conserva- 
tion law as follows from (|5.5|) , it may also admit new torsion-dependent integrals of motion 
0. To illustrate our statement, consider two-dimensional motion in the constant metric 
and torsion fields, S^ vc = tT w , where T va = —T av is the generator of SO(2), T i2 = 1, 



but this symmetry does not lead to the conservation of the momentum components as 



and d^ u = 0, = 5^. The Lagrangian L = vjl exhibits the translational symmetry, 



one might see from (|5.5 ). Nonetheless, we may solve Eq. ( |5.6| ) relative and find new 



integrals instead of the Noether's integrals. The solution is $ = 0,lu^ = [exp((pT)]^a v 
3 This conclusion seems to be opposite to the one made in |15|. 
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where a M is an arbitrary constant vector and p = 2<5 / ^7 M g M . Since a M is arbitrary, we have 
two independent integrals of motion 

l\ = V\ cos p — v 2 sin p , I 2 = V\ sin p + v 2 cos p . (5.7) 

The integrals ( p.7[ ) comprise two independent first integrals for the two-dimensional au- 
toparallel motion in the homogeneous metric and torsion fields. We remark also that 
I\ + If = IE where E = L is the energy. 

Instead of modifying the transformation law of the Lagrangian, one can modify the 
transformation law of its arguments, the generalized coordinates and velocities. Set 

d w q» = u» , duv" = d v u" + 25" w w v / . (5.8) 

If the Lagrangian is invariant under these transformations up to a total time derivative (see 



(|5.2| )), then (|5.3|) is an integral of motion. The attention should be paid to the fact that 
in contrast to the conventional formulation of Noether's theorem [HJ, the transformation 
law (|5.8| ) determines no smooth deformation of the path q fl (s) on a manifold and, in 
this sense, is not induced by any diffeomorphism on M. In fact, the relations ( |5.8| ) is 
a postulate that specifies the transformation law of two independent coordinates of the 
tangent bundle. One first calculates d w L(v,q), then one sets t> M = if and looks for such 
tjjv = uj^(v, q) that ( |5.2|) holds for some $. This procedure, though being somewhat 
unusual, may turn out to be useful in seeking the integrals of motion for the modified 
Euler-Lagrange equations ( |4.11| ). 

As an example, consider the autoparallel motion in the teleparallel space (zero Rie- 
mann-Cartan curvature and non-zero torsion). Such spaces are used to describe a crystal 



with topological defects [17|. In the anholonomic system discussed in Section 2, we set 
fuv = 0' then the curvature vanishes, while the torsion does not. We reduce either of the 
actions ( |1.3| ) and ( |1.4| ) on the constraint surface and obtain L = L(v,q). Now it is not 
hard to verify that the transformations (|5.8|) with = (e 11 , a), a 1 arbitrary constants, 



leaves the Lagrangian invariant ($ = 0). Since a 1 are arbitrary, the quantities 

r = 4( ? K (5-9) 

are the integrals of motion. Thus, the autoparallel motion in the teleparallel spaces can be 
characterized by a simple property: The velocity components taken in the non-coordinate 
basis [ti = efd^ in T q M) that is transported parallel {D^e % v = 0) are the integrals of 
motion (P = 0). This is, in general, not the case for non-teleparallel spaces. 



6 A Lagrangian formalism for autoparallels 

As has been pointed out in the introductionary remarks, there might exist a non-local 
action whose extrema (relative to the conventional variational principle) would contain 
autoparallels. In fact, there are infinitely many such actions. We shall give one of the 
possible actions. It has a few additional properties that seems to us useful for the canon- 
ical quantization of the autoparallel motion. We require that the sought action should 
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coincide with a local action whose extrema are geodesies when the torsion is zero, and 
the non-locality can be removed by adding new degrees of freedom coupled to the original 
variables. Thus, the autoparallel motion can be modeled by a holonomic dynamics with 
some auxiliary degrees of freedom which admits both a Hamiltonian formalism and the 
canonical quantization. 

We skip the details of our derivation and just give the answer. Let S g be a local action 
whose extrema are geodesies, i.e., we set 

H = 9^ . (6.1) 
Let us introduce an integral operator A by the relation 

J ds'A^s, s') D v v v {s') = 9fl J(s - s") . (6.2) 

Here and below the integration is extended from the initial to final moment of time. 
Consider the non-local action 

S = S g + Jds (D v v» - D v v>*) A^D v v u = S g + J dsK% x v»v x A ua D v v a . (6.3) 

Making use of the properties ( |6.1| ) and ( |6.2[ ), one can convince oneself that from D v v^ = 
follows 6 q ^S — 0. So the autoparallels are extrema of the action (|6.3| ). If the torsion is 
zero, the non-local term in ( |6.3|) vanishes and S = S g . 

The non-locality of the action (|6.3| ) still prevents us from developing a Hamiltonian 
formalism. We need a local action. Fortunately, the action (|6.3p can be regarded as an 
effective action for physical degrees of freedom in a larger dynamical system which is 
described by a local action. Let us extend the original configuration space of the system 
by auxiliary variables y^ whose dynamics is determined by the equation 

y" = g^A va D v v° . (6.4) 

Consider the action 

S = S g + Jds [^«Y + A M (sTtlnif - D v v»)] . (6.5) 



The operator Vt is to be chosen so that the action ( |6.5|) turns into the action fl6.3|) when 
it is reduced on the solution of the equation of motion for the variable A M . 
We have 

5S r 

9M = J ds'to^s, s')y°(s') - 9iiV D v v v {s) = . (6.6) 

From Eq. fl6.2|) follows that the kernel of the operator Cl can be taken in the form 

tt^(s, s") = g^(s) Sga ,„j D v v u (s) . (6.7) 
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With this choice Eq. ( |6.6| ) becomes an ordinary linear differential equation of second 
order for y M : 

r + + r" w ) v°y u + y v d v Y\ x v°v x = D v v» . (6.8) 

Its solution is a sum of a general solution to the homogeneous equation Cl fJiU y u = and 
a special solution y^ = g^ v K va D v v a . To recover the action (|6.3|) , we have to supplement 
Eq. ( |6.8| ) by the initial conditions such that the homogeneous equation has only a trivial 
solution. This is always possible thanks to the linearity of the equation. In particular, 
one can take y^ = y^ = at the initial moment of time. 

The local action (|6.5| ) linearly depends on the second-order derivatives y^ and q**. 
Assuming zero boundary conditions for the variable A M , we may remove the second-order 
time derivatives by integrating by parts, thus producing the final local Lagrangian that 
depends only on the coordinates and velocities in the extended configuration space and 
involves no higher-order time derivatives. It is important to observe that the Hessian 
of the Lagrangian is not degenerate, that is, the system exhibits no constraints. The 
dynamics of the auxiliary degrees of freedom y M and A M has been chosen so that when 
the torsion is present, the coupling between them and the original variables g M causes the 
deviation of the trajectory g M (s) from the geodesic, making it the autoparallel. 

In the path integral formulation of quantum mechanics, the auxiliary degrees of free- 
dom A M and y^ must be integrated out to obtain an effective path integral for the original 
system. This resembles the Feynman- Vernon approach to quantum dissipative systems 



181 , where the non-potential friction force is generated by a special coupling of the oscil- 
lator bath to the system. In our case, the non-potential (non-Lagrangian) torsion force is 
modeled by a special coupling to the A- and ^/-degrees of freedom. Extending this anal- 
ogy further, one may expect that quantum mechanics in spaces with torsion, that favors 
autoparallels in the classical limit, can not, in general, be described by the wave function 
formalism, rather only the density matrix can be constructed like for dissipative systems. 
The nice property which could be expected is that in the semiclassical approximation 
the transition amplitude generated by the effective action ( |6.3| ) is given by the geodesic 
action S g taken on the autoparallel, the non-local term of (|6.3|) will contribute only to 
quantum fluctuations because D v v M = 0. The canonical quantization of the model will be 
considered elsewhere. 



7 Conclusions 

We have analyzed the autoparallel motion from the point of view of analytical mechanics 
and succeeded to represent it as a very special anholonomic constrained system. Invoking 
the variational principles for anholonomic dynamical systems, we have established the 
covariant variational principle for the autoparallels. We have also analyzed a modification 
of Noether's theorem due to the torsion force. Finally, we have found a possible local action 
whose extrema determined by the conventional Hamilton variational principle contain the 
autoparallels for some degrees of freedom. The model can be canonically quantized by 
means of the Dirac formalism [Q]. 
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